Domain- Wall Waves (2D Magnons) in Superconducting Ferromagnets 
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Propagation of the magnetization waves along domain walls (2D magnons) in a superconducting 
ferromagnet has been studied theoretically. The magnetostatic fields (long-range dipole-dipole inter- 
action) have a crucial effect on the spectrum of 2D magnons. But this effect is essentially affected 
by the superconducting Meissner currents, which screen the magnetostatic fields and modify the 
long-wavelength spectrum from square-root to linear. The excitation of the domain wall waves by 
an electromagnetic wave incident on a superconducting-ferromagnet sample has been considered. 
This suggests using measurements of the surface impedance for studying the domain wall waves, 
and eventually for effective probing of superconductivity-ferromagnetism coexistence. 

PACS numbers: 74.25.Ha, 74.90.+n, 75.60.-d 



o 
o 

Of 
C3 



-a 
c 

o 
o 



> 

OS 

as 

SO 

o 

so 
o 

-I— > 



I 

c 

o 
o 



X 



Coexistence of superconductivity and ferromagnetism, 
which results in a number of unusual phenomena, has 
already been studied about 50 years A revival of 
research in this area was stimulated by experimental 
observation of superconductivity-ferromagnetism coexis- 
tence in high-T c superconductors I2U3J and various un- 
conventional superconductors ^OjlTQ- An essential 
obstacle for experimental detection of superconductivity- 
ferromagnetism coexistence is screening of internal mag- 
netic fields generated by the ferromagnetic order parame- 
ter (spontaneous magnetization) by superconducting cur- 
rents. A possible way to overcome this obstacle is to in- 
vestigate spin waves, which are a direct evidence of the 
presence of spontaneous magnetization 0, Q . They can 
be excited by the electromagnetic (EM) wave incident on 
a superconducting ferromagnet (SCFM). 

In Ref. y spin waves in SCFMs have been studied for 
the Meissner state with the uniform spontaneous mag- 
netization, i.e., for a single-domain sample. Though at 
the equilibrium superconductivity suppresses usual fer- 
romagnetic domain structures with periods determined 
by demagnetization factors and sizes of samples |lf| , the 
interplay with superconductivity can lead to formation of 
intrinsic domains of the order of the London penetration 
depth or less (see the latest discussion in Ref. HaV 
Moreover, a domain wall (DW) is a topologically stable 
defect, and the presence (absence) of DWs can be deter- 
mined by the sample prehistory but not by the conditions 
of the equilibrium. This justifies an interest to studying 
spin waves propagating along DWs (2D magnons). 

In normal ferromagnet (FM) the DW dynamics and 
DW waves have already been studied a few decades and 
are important for various applications [lflj . The crucial 
feature of 2D magnons is that the magnetostatic effects 
(dipole-dipole interaction) are much more important for 
them than for 3D bulk magnons. On the other hand one 
may expect that SC should influence these effects first 
of all since it prohibits penetration of the magnetostatic 
fields deep into the bulk of domains. 

In the present Letter we investigate theoretically DW 



waves in a SCFM. The analysis confirms expectation of 
strong effect of superconductivity on DW waves: Meiss- 
ner screening of magnetostatic fields in domains modifies 
the square-root 2D-magnon spectrum revealed for nor- 
mal FM 0| to a linear sound-like spectrum. We show 
how the 2D magnons can exhibit themselves in the linear 
response to the AC magnetic field (surface impedance). 
As well as in the case of bulk 3D spin waves, detection 
of 2D magnons, being interesting itself, would provide a 
direct probe of superconductivity-ferromagnetism coex- 
istence. The spectrum of 2D magnons is sensitive to bulk 
and surface pining of DWs, which can lead to a gap in 
the spectrum and a peak in the EM wave absorption. 

Let us consider a 180° DW in a FM with magnetic 
anisotropy of the easy-axis type [Fig.^a)]. The z axis is 
the easy axis and the DW is parallel to the xz plane and 
separates two domains with spontaneous magnetization 
M along the +z (y < 0) and — z (y > 0) directions. 
Neglecting the magnetostatic fields, the structure of the 
DW is well known [llQ- Rotation of the magnetization 
M = M (cos<f) sin 6*, sin sin 9, cos 9) inside the DW is 
described in polar and azimuthal angles 9 and <j> as 



(9n = 2 tan 



-1 e y/A 



(1) 



Here A = (A/K^^Mq is the DW width, A is the con- 
stant, which determines the energy of the non- uniform 
exchange v4VjM • V^M, and the constant K determines 
the easy-axis anisotropy energy if sin 2 9. The azimuthal 
angle <j) is an arbitrary constant as far as the magneto- 
static (dipole) fields are neglected. In normal FMs cur- 
rents are absent at the equilibrium , and the magnetic 
field H = B — 47rM is curl- free, i.e., is a potential field 
with sources called "magnetic charges" p m = — V • M: 
V • H = V ■ B - 4ttV ■ M = 47rp m . The energy of the 
magnetostatic fields lifts degeneracy with respect to the 
angle <j>: the ground state corresponds to the DW of the 
Bloch type with <j> = 0. This means that the magnetiza- 
tion rotates in the DW plane, and the magnetic charges 
do not appear (V • M = 0). This structure, which was 
obtained for normal FMs, remains valid at scales of or- 
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FIG. 1: (Color online) (a) The 180° Bloch domain wall. The 
magnetizations in the domains are parallel and antiparallel to 
the z axis. Inside the DW the magnetization rotates in the xz 
plane, (b) Excitation of DW waves by an incident EM wave. 



der A also for SCFMs, as far as the London penetration 
depth A essentially exceeds the DW thickness A. How- 
ever, the difference between a normal FM and a SCFM 
is important at distances larger than A: while in nor- 
mal FMs the magnetic field H — B — 47rM vanishes and 
the magnetic induction B = 4-7rM is constant inside do- 
mains, in SCFMs the magnetic induction B is confined 
in the Meisssner layers of width A [Hj : 



B z = ±4:nM e ±y/x 



(2) 



where the upper and the lower signs correspond to y < 
and y > respectively. Thus the Meissner currents 
j x = —(c/A^dBz/dy screen out the main bulk of do- 
mains from the magnetic induction. 

We shall use the DW dynamics developed for FMs with 
high quality factor a = K/2itMq 0|. In this limit one 
may assume that the structure of the Landau-Lifshitz 
DW, which is given by Eq. (J), is not affected essentially 
by dynamical processes, and the DW dynamics can be 
described in the terms of the pair of canonically con- 
jugated variables: the azimuthal angle (f>(x, z;t), which 
determines the plane of magnetization rotation, and the 
displacement of the DW r\ — r)(x,z;t) along the y axis: 
9(x,y,z;t) — 9o(y — rj(x, z; t)). Neglecting dissipation, 



the equations of DW motion (Slonczewski's equations 
[Tq) are the Hamilton equations: 

SH „ , 7 SH 
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2M 



d t <t> = 



(3) 



2M Srj ' 

where 6 /Sip and S/5rj are functional derivatives and the 
Hamiltonian H = J adx dz is determined by the surface 
energy density a of the DW. The Hamilton equations © 
point out a direct analogy of (2Mo/j)(f> with a canonical 
momentum conjugate to "coordinate" rj. But it is worth- 
while to mention another analogy. The displacement rj 
leads to the change 2Mqt\ of the z-componcnt of the to- 
tal magnetization per unit area of the DW. Since <p is 
the angle of spin rotation around the z axis, Eqs. J3J can 
be also considered as Hamilton equations for the pair of 
conjugate variables "angular momentum-angle" . 
The DW surface energy is 

a{(j),rj) = 47rM 2 Aa [(Vr/) 2 + A 2 (V0) 2 ] 
+i- J dy(h 2 + A 2 [Vxh] 2 ) 

+ ^ J dy(H<® • h 2 + A 2 [V x H<°>] • [V x h 2 ]) . (4) 

The terms oc a (the first square brackets) present the 
contributions of the exchange energy and the anisotropy 
energy, whereas the integral terms are the energies of 
superconducting screening currents and of the magneto- 
static fields H = H(°)+h+li2. Here H^ ^ is the magnetic 
field in the ground state, and h and h 2 are the dynamic 
corrections to the field of the first and the second order in 
the wave amplitude respectively. Necessity to take into 
account the second-order corrections is a peculiar feature 
of the DW dynamics in SCFM: in normal FM H<°) = 
and these corrections vanish. 

The first-order contribution to the magnetic induction 
b = B — B'°) = h + 4nm can be found from the gener- 
alized London equation: 



A b + A" 2 b = 4?rV x V x m. 



(5) 



The dynamic components of the magnetization m = M 

arc defined through the conjugated variables 77 and 
6 as follows: 



m = Mr 







-7yd[sin 9 (y)]/dy\ 

0sin6» o (y) 
-(rj/A)sm 2 e (y)J 



(0) 



Let us consider the plane wave with wave vector k = 
(0, 0, k) parallel to the z axis: 0, 77, b oc exp(— iuit + ikz). 
Then the London equation (0) becomes a set of ordinary 
differential equations with the only coordinate y, and its 
solution is: 



47T 

b = — 

k 



+e 



-ky 



dy'V x V x m(y')e 
dy'V x V x m(y') 



-ky' 



o k v' 



(7) 
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Here fc 2 = fc 2 + A~ 2 . 

Now one can substitute the magnetic field h = b— 47rm 
into the first magnetostatic integral in Eq. J3J and 
perform integration. The contribution of the second- 
order corrections h.2 [the second magnetostatic integral 
in Eq. Q]can be determined from general arguments. 
The second-order terms contain the Fourier components 
of rj 2 with zero wave vector and with the wave vector 
Ik (the second harmonic). Only the first component, 
which is z-independent, can interfere with the zero-order 
terms. This yields the fc-independent contribution to 
the surface energy, which should provide translational 
invariance: the surface energy should not depend on the 
spatially independent DW displacement. Altogether this 
yields the following surface energy density in variables ry 
and 6: 



1 - 



a = 4nM 2 [afc 2 A(77 2 + c/> 2 A 2 ) 



Afc 2 



l A- 



(8) 



where the second-order terms are presented by the nega- 
tive term — r/ 2 /A. All terms, which are not proportional 
to the quality factor a, originate from the magnetostatic 
energy and the SC currents. The magnetostatic term 
oc (f) 2 plays a role of the "kinetic" energy since the angle 
4> is the momentum conjugate to the displacement 77. In 
the limit fc — > it is determined by the magnetic fields 
inside the DW and can be expressed via the Doring mass 
mo = l/7 2 27rA ^jj. The fc-dependent contribution to 
the kinetic energy represents the energy of the magneto- 
static fields outside the DW in the layer of width ~ l/|fc|- 
Altogether the magnetostatic terms oc 4> 2 are connected 
with the double layer of magnetic charges, which appears 
when the DW rotates from the Bloch orientation = 0. 
On the other hand, the magnetostatic term oc rf 2 is con- 
nected with the magnetic charges, which appear due to 
rotation of the DW with respect to the direction of the 
magnetization in domains. 

Using Eq. © in Slonczewski's equations 10 trans- 
formed to the Fourier presentation one obtains the spec- 
trum of plane waves propagating along the DW: 



afc 2 A 2 + j (Vfc 2 A 2 + 1 - l) 
fc 2 AA 



1 



*fc 2 A 2 - 



v / fc 2 A 2 



1 



(9) 



where u>m = ^irjM. 

In the limit A —> 00 (no superconductivity) , this trans- 
forms to the spectrum of 2D magnons in normal FMs: 

to 2 = u 2 M (ak 2 A 2 + |fc|A)(l + afc 2 A 2 - |fc|A). (10) 

All nonanalytic terms oc |fc| in this expression are of mag- 
netostatic origin and are related with penetration of mag- 
netostatic fields deep into domains on scales of the order 



of the wavelength 2ir/k. They lead to the nonanalytic 
wave spectrum uj oc \f~k in the long- wavelength limit, 
which was known before [T^j ] . Another nonanalytic term, 
which appears in the second multiplier with the nega- 
tive sign, has not been considered so far as far as we are 
aware. 

Returning back to SCFM one can see that the Meissner 
screening eliminates non-analytic features of the spec- 
trum and allows expansion in fc 2 . Thus in the long- 
wavelength limit the spectrum is sound-like: 

uj = c s k, (11) 



where c s = wi(\/A(aA + A/2) is the 2D spin- wave ve- 
locity. 

In order to find the effect of DW waves on surface 
impedance one should solve a boundary problem. In gen- 
eral the wave spectrum Eq. © leads to differential equa- 
tions of high order in the configurational space (bear- 
ing in mind the correspondence fc — > —id/dz). But the 
sound-like spectrum Eq. i|llfl reduces Slonczewski's equa- 
tions in the configurational space to two ordinary differ- 
ential equations of the first order: 



dtV = wmA^ - jAh ext , d t 4> ■ 



d 2 r\ 



uj m A dz 2 



(12) 



In the first equation we added the interaction with the 
external magnetic field. This field arises from the EM 
wave of frequency uj incident on the sample surface and 
linearly polarized along y-axis [Fig.^b)]. The EM wave 
penetrates through the surface of SCFM at the distance 
~ A: h ext = h$e z l x . Here ho is the AC magnetic field at 
the sample surface z — 0. 

Next is to formulate boundary conditions for these 
equations. One boundary condition is imposed at z — > 
—00. We assume that the DW plane wave, which is ex- 
cited near the sample surface z — is the only propa- 
gating wave in the bulk and there is no reflected wave 
coming from z — — 00. The second boundary condition 
is imposed at the sample border z = 0. Assuming an 
ideal surface without any surface force the balance of the 
"momentum" (f> requires that the DW remains normal to 
the sample border: dr]/dz = 0. But the magnetic field 
h ext generated by the incident EM wave , though being 
a bulk force, is confined to the Meissner layer of width 
A, which is much smaller than the wavelength. Then one 
may consider its integral effect as that of a surface force. 
This leads to modification of the boundary condition for 
propagating wave, which must be written as 



,dr) 
'dz 



z=0 



— iuj^h^AX = 



(13) 



This boundary condition determines the amplitude of the 
wave propagating from the sample border to z — — 00: 



V 



-7/ioAAexp (—iu>t — ikz) 



(14) 
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The wave is accompanied by the energy flux, which is 
determined from the energy balance dtcr + dS/dz = 0: 



S=- 



da 



d{drj/dz) 



16 c. 



2,2 
"0 



(15) 



The energy brought away by the DW plane wave is sup- 
plied by the wave incident on the surface: at the reflection 
of the wave from the sample surface some part of energy 
is absorbed though no dissipation mechanism is present 
explicitly in our model: it is assumed that dissipation 
occurs deep inside the sample where the spin wave even- 
tually dissipates and does not return as a reflected wave. 
Energy absorption at reflection of the EM wave from the 
sample surface is characterized by real part of the surface 
impedance, which is defined as a ratio of the tangential 
components of electric field e x and magnetic field ho at 
the surface: £ = e x /ho\ z =o- Namely, absorption per unit 
area per second is given by (ch^ /8n)Re(. Equating this 
to the average energy flux per unit area given by |if?|niy 
one obtains that 



I CI 2 2 

ReC = Snnw i-i = ^L. nw AX 2 — 
cni 2 cc* 



(16) 



Here nw is the linear density of DWs. However this 
dependence is essentially modified by bulk and surface 
pinning of DW. Bulk pinning lifts translational invariance 
and leads to the gap in the magnon spectrum, which 
instead of Eq. l(TT)) becomes 



2 i 2; 2 

UJ p + c s k 



(17) 



Surface pinning modifies the boundary condition 
Eq. i|13|) . which now is 



oz 



iwyhoAX = 



(18) 



z=0 



where k s is the inverse length characterizing intensity of 
surface pinning. Repeating the derivation of ReC in the 
presence of pinning one obtains that 



ReC = — n w AA 2 — —= 
I cc s (ur 



aA/w 2 - ul 



(19) 



Thus in the presence of bulk pinning absorption of the 
EM wave starts from the threshold frequency u> p and 
above the threshold grows as y/ui — oj p . If at the same 
time surface pinning is weak enough there is an absorp- 
tion peak at u> — lo p ~ c s k s . The peak can correspond 
to frequencies much lower than the threshold frequency 
for excitation of bulk spin waves (frequency <xu>m of the 
ferromagnetic resonance) These special features of 
the absorption spectrum can be useful for experimental 
detection of DW waves in SCFMs. 

Our analysis was based on the classical Landau-Lifshitz 
dynamics of a spin FM with a single- valued spontaneous 



magnetization M. Meanwhile in p-wave superconductors 
superconductivity coexists with orbital ferromagnetism 
when spontaneous magnetic moment is not well-defined 
and one should develop the magnetic dynamics using the 
concept of magnetization currents |lg. But as well as 
in the case of bulk magnons considered in Ref. 0, one 
may expect that this would lead only to redefinition of 
the 2D magnon parameters without changing the general 
picture of their propagation and excitation. We hope to 
address this problem elsewhere. 

In summary, the low-frequency dynamic of the DW in 
SCFM has been studied and the spectrum of DW waves 
(2D magnons) has been found. The magnetostatic effects 
(long-range dipole-dipole interaction) have a crucial in- 
fluence on the long-wavelength spectrum, while the su- 
perconducting Meissner currents modify these effects es- 
sentially. They suppress the effect of long range dipole 
interaction eliminating the non-analytical features of the 
2D-magnon spectrum revealed for normal FMs. Our 
analysis has demonstrated that the response of a SCFM 
sample to the EM irradiation (surface impedance) pro- 
vides information on 2D magnons in DWs. Bulk and 
surface pinning of DWs introduces a gap in the 2D- 
magnon spectrum and a characteristic peak in the ab- 
sorption of the incident EM wave (real part of the sur- 
face impedance). In conclusion, studying of DW waves 
promises to be an effective probe of superconductivity- 
ferromagnetism coexistence. 
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